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Abstract. The Lagrangian proposed by York et al. and the covariant first 
order Lagrangian for General Relativity are introduced to deal with the (vac- 
uum) gravitational field on a reference background. The two Lagrangians are 
compared and we show that the first one can be obtained from the latter under 
suitable hypotheses. The induced variational principles are also compared and 
discussed. A conditioned correspondence among Nother conserved quantities, 
quasilocal energy and the standard Hamiltonian obtained by 3 + 1 decompo- 
sition is also established. As a result, it turns out that the covariant first order 
Lagrangian is better suited whenever a reference background field has to be 
taken into account, as it is commonly accepted when dealing with conserved 
quantities in non-asymptotically flat spacetimes. As a further advantage of 
the use of a covariant first order Lagrangian, we show that all the quantities 
computed are manifestly covariant, as it is appropriate in General Relativity. 

1. Introduction 

Many approaches to variational principles, conserved quantities and related 
topics can be found in the current literature about General Relativity (see, e.g., 
[1], [2], [3], [4], [5], [6], [7], [8], [9], [10], [11], [12], [13] and references 
quoted therein). We shall hereafter compare two of them, both dealing with 
a dynamical metric g and a reference background metric g over a spacetime 
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manifold M of arbitrary dimension n , with n > 3 . A large part of the paper 
will deal with the explicit case n = 4 , in view of applications to standard 
General Relativity. We shall deal with the vacuum case; the generalization to 
the case of presence of matter fields is straightforward. 

The first variational principle we shall deal with is based on the covariant 
first order action functional for General Relativity (see [12], [14], [15], [16]). 
The second one ( [3], [7]) is based on the action functional due to York et al., 
which is defined with the aim of dealing with the fixing of the (n — 1) -metric 
induced on the boundary dD of any region D of spacetime M . 

The covariant first order action was introduced to set General Relativity 
in a standard covariant first order variational framework. As is well known, 
in fact, the Hilbert action functional L = (1 /2k) ^fg Rds is second order in 
the metric field, so that field equations are expected to be of the fourth order. 
Einstein field equations are second order equations instead, as if the action 
were first order only. This is due to the well known fact that locally second 
derivatives of the metric field appearing in the scalar curvature may be hidden 
under a divergence, thus not appearing in field equations (a fact which was clear 
to Einstein from the very beginning; see, e.g. [17]). Of course, however, this 
cannot be done in general in a global and covariant way; that is why General 
Relativity is usually considered as a second order field theory or, whenever it 
is treated on a first order basis, something is lost (e.g. covariance or boundary 
terms; see [4], [5]). 

The covariant first order action functional is the following: 

A ° [9r9]=i LJ D v ° Rds ~ i kJ dD v^"x„ ds «-^J D ^ Rds 

where k is a constant ( k = 8ttG/c 4 in General Relativity with dim(M) = 4 ), 
^fg is the square root of the absolute value of the determinant of the dynamical 

metric g , ^/g = ^/|det^| the analogous quantity for the background metric 
g and ds and ds a = <9 a jds are the standard local bases for n -forms and 
(n — 1) -forms over M , respectively. We systematically denote by a bar the 
quantities referred to the background, i.e. here and hereafter we shall use the 
following notation 
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Let us also introduce the following relative quantities: 

'iaP — L a/3 L a/3 > W a/3 ~ U ci(3 U ci(3 



The action functional (1.1) is associated to the so-called covariant first 
order Lagrangian 

L=±(ygR-d Q (V&T - R) ds (1.2.a) 

or equivalent ly 

From the first expression (1.2. a) , it can be easily seen that the fields g and g 
do not interact and they both obey vacuum Einstein field equations (provided 
that suitable boundary conditions are satisfied, as we shall discuss shortly after). 
From the second expression (1.2.6) , the Lagrangian is however recognised to 
be first order in g and second order in g . Being both q^ and tensors, 

L is a covariant Lagrangian. It is, of course, the truly covariant counterpart 
of the so-called Hilbert-Palatini first order Lagrangian (see [18]), to which it 
reduces by suitable non-covariant cancellations and background fixings. We 
stress that in the variational principle induced by (1.1) the dynamical metric 
g is endowed with a direct physical meaning, while the reference background 
metric g is, at least for the moment, introduced to provide covariance and as a 
reference value for conserved quantities, as discussed below. Notice that if the 
action is computed for g = g , then it identically vanishes. 

We remark that the Hilbert-Einstein Lagrangian L = (I /2k) ^fg Rds is 
a second order Lagrangian and thus it would a priori define a variational prin- 
ciple in which the dynamical metric g together with its first derivatives are 
kept fixed on the boundary. In fact, in standard variational calculus, for a co- 
variant field theory described by a Lagrangian of order k , one keeps fixed the 
fields on the boundary together with their derivatives up to order k — 1 (see 
Appendix B). However, the degeneracy of the Hilbert-Einstein action func- 
tional has been globally removed by introducing a reference background metric 
g which provides a covariant way to cancel out the second order derivatives of 
the dynamical metric g , as shown in (1.2.6) above. As we already remarked, 
the Lagrangian (1.2.6) is first order in g and second order in g so that it de- 
fines a standard variational principle in which the dynamical metric g is kept 
fixed on the boundary, while the reference background metric g is kept fixed 
together with its first derivatives on the boundary. Of course, nothing more 
can be done without breaking down covariance, i.e. without choosing an ADM 
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foliation, or without requiring suitable additional matching conditions between 
the two metrics. 

On the other hand, York's action was originally built (see [3], [7], [19]) 
to provide a variational principle suited to deal with the boundary conditions 
which are specified by keeping the induced metric on dD fixed. Boundary 
terms are added to the standard Hilbert-Einstein action (1/2k) f D R^/gds. 
They are exactly needed so that the boundary contribution to the variation of 
the action vanishes when the induced (n — 1) -metric is kept fixed on dD . 

The York's action functional is adapted to an ADM foliation induced by 
dragging a spacelike hypersurface E along a timelike vector field ( . 

An ADM foliation of a spacetime re- 
gion D is obtained and it is parame- 
trized by the affine parameter t of £ 
(see Fig. 1 ). The region D is thence 
topologically the product of E times 
a real line interval [to>£i]- Let the 
generic leaf be E t and B t its bound- 
ary, which is obviously (n — 2) -dimen- 
sional. Let us then denote by u the 
future directed timelike unit normal to 
the leaf E t . The boundary dD is 
formed by the union of all (n — 2) - 
boundaries B t , which will be denoted 
by B , together with the initial and fi- 
nal leaves E to and E tl of the sand- 
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wich, which will be called lids. Let 
us denote by n the outward spacelike 
unit normal to B , by 7^ the metric 
induced by g on B , by h a b the metric induced on E t and by ab the metric 
induced on the (n — 2) -boundary B t . Here and in the sequel indices run in 
the following ranges: Greek indices from to n — 1 , lower case Roman letters 
a, 6, . . . from 1 to n — 1 , lower case middle Roman letters i,j, . . . take the 
values 0,2, . . .n — 1 and upper case Roman letters A, B . . . range from 2 to 
n — 1 (see also Appendix A for the notation). 

We shall finally denote by 0y the extrinsic curvature of B in the space- 
time M , by K a b the extrinsic curvature of E t in M and by /Cab the extrinsic 
curvature of B t in E t (see Appendix A). We shall denote by = 7 lJ , 
K = h ab K a b and /C = a AB /Cab the traces of the extrinsic curvatures of B , 
Et and Bt , respectively. It is assumed that the dynamical metric g and the 
reference background g induce the same metric 7^ on B . Furthermore, let 
us also assume that the hypersurfaces B and E t intersect orthogonally (or 
equivalently that n and u are orthogonal on any B t , i.e. -u M n^\s = ). We 
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stress that here we are considering a region D in which the timelike vector field 
£ has no fixed points, so that the hypersurfaces S t do not intersect each other 
and span the whole region D . 

According to this notation, York's action in presence of a background g 
may be written as 

I D [g,g]=I D [g}-I D [g} (1.3) 
where the functional Id [g] is defined by: 

Id[9\=^- [ ^gRds + - / H Ku a ^gds a -- [ Qn a ^/gds a (1.4) 

*K J D K JH to K JB 

and Id [g] is the same functional calculated for the background g . Notice that 
y/gu a ds Q = \/hd s x and y^n a ds a = d 3 x are the volume elements on E t 
and B , respectively. In the functional (1.4) we also set the convenient notation 

•^£ to "fctx "^t 

In the current literature, there exists a whole family of action functionals 
similar to (1.3) each adapted to the particular problem under consideration. In 
general, one can add to the functional (1.4) an arbitrary functional depending 
on the data fixed on the boundary, i.e. depending on the boundary metric 
(see [3]). This arbitrariness does not affect the equations of motion since the 
boundary metric is kept fixed in the variational principle. The choice (1.3) is 
motivated by the requirement iu[<7, g] = , i.e. the requirement that the action 
functional vanishes when computed with g = g . The same property is satisfied 
by (1.1) • 

Moreover, under additional hypotheses, the integral on the lids T, to and 
E tl in (1.4) is usually discarded since the lids are either ignored or identified 
by periodic boundary conditions (see [20], [21]). 

For these reasons we are forced to first review the literature and choose 
homogeneous notation to compare with the covariant first order Lagrangian. 

We are not aware of a theoretical comparison between these two varia- 
tional principles in the current literature. A definition of energy based on the 
action functional (1.3) has been proposed by Brown and York. It is called 
quasilocal energy and it is often quoted because it reproduces the ADM mass 
in the asymptotically flat case, though it has been recognised to be also suit- 
able for more general boundary conditions (e.g. in asymptotically anti-de-Sitter 
(see [22]) or asymptotically locally flat cases; see [23]). The main advantage 
of Brown- York method is that it allows to define the energy within a finite 
region as well as the total one. We remark that the covariant first order La- 
grangian (1.1) was originally proposed for exactly the same reasons (see [14], 
[15]). We also remark that in both cases it has been recognized that absolute 
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conserved quantities have a meaning just in particular cases, while for general 
boundary conditions just conserved quantities relative to a reference background 
are meaningful. Furthermore, even "absolute" conserved quantities should be 
interpreted as conserved quantities relative to some canonical background (e.g. 
flat Minkowski space in asymptotically flat spaces). We also remark that refer- 
ence backgrounds are particularly relevant in General Relativity, as well as in 
other non-linear field theories, since whenever fields are endowed with a vec- 
tor space structure then a canonical choice for the reference exists, namely the 
zero section. If the configuration bundle (see Appendix B for a short geometric 
insight into variational calculus and field theory) is not a vector bundle, as it 
happens in General Relativity, as well as e.g. in Yang-Mills theories, there is no 
canonical choice for the vacuum state. The vacuum state has thence to be ar- 
bitrarily fixed. When doing that it sounds physically reasonable to require also 
the background to be a solution of field equations, so that the relative mass can 
be interpreted as the energy "spent to go" from the background solution to the 
dynamical one (and analogously for other currents). Furthermore, it is essential 
that the choice of the reference background does not effect the evolution of the 
dynamical fields, i.e. they have to be decoupled. In view of these considerations 
both action functionals (1.1) and (1.3) incorporate the background from the 
very beginning. Thus we believe that the relationship between the two methods 
is worth investigating. 

In this paper we shall also discuss the very notion of conserved quantity. On 
one hand, in fact, Nother theorem provides currents £ which are covariantly 
conserved, i.e. dS = identically or on shell (i.e., along solutions), meaning 
that their integral on the boundary dD of any n -region D in spacetime 
M vanishes or, equivalently, that the conserved quantity obeys a continuity 
equation. On the other hand, physicists are often interested in quantities which 
are conserved in time, meaning that, once an ADM foliation of a region of 
spacetime has been chosen, such a quantity Q may be computed by integration 
on each leaf and it turns out not to depend on the particular leaf labelled by the 
time t . Clearly, ADM foliations are far not unique and different foliations of 
the same region D correspond to different ways of defining time. Furthermore, 
such a quantity Q may be conserved in the time defined by an ADM foliation 
without being conserved with respect to other foliations (see Appendix C). From 
a theoretical General Relativity viewpoint, quantities conserved in time are not 
(manifestly) covariant in nature. They are, in fact, conserved with respect 
to a special parameter, while, at a fundamental level, the principle of general 
covariance forbids, at least in principle, the selection of a preferred time. 

Nevertheless, quantities conserved in time may be interesting to be inves- 
tigated. In our perspective, in fact, they can be obtained from covariantly 
conserved quantities. To be more precise, we can consider a variational prin- 
ciple, an infinitesimal generator £ of Lagrangian symmetries and a solution of 
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field equations. Then we can compute covariantly conserved currents £[£] by 
Nother theorem. Let us then fix a spacelike (n — 1) -region E and integrate 
the Nother current on it to define a conserved quantity . Any timelike 
vector field ( allows then to evolve the region E along its flow, parametrized 
by its affine parameter t . Under this viewpoint, the question arises whether 
there exists a vector field £ (possibly depending on the region E ) such that the 
covariantly conserved quantity generated by £ is also conserved in the "time" 
induced by ( . At a first glance, if we have dD = E tl — E to + B (see Fig. 
1), conservation in time is equivalent to require that the integral of the Nother 
current £[£] on B vanishes (for any time interval [£o>£i] )• In fact, we have 
the covariant conservation law d£[£] = , so that 



o = / = [ m = [ m - [ m + [ m 

JD JdD JY, tl Jn t0 J® 

[ m - [ m = - [ m 



(1.5) 



Thence the conserved quantity / s £[£] computed on a leaf does not depend 
on the particular leaf if and only if J B £[£] = (for any time interval [to, t±] ). 
Physically speaking, this amount to require that the flow of the current £[£] 
through B is vanishing. Clearly, different ADM foliations may evolve E in dif- 
ferent ways. In general, just few of them will lead to time-conserved quantities. 
The vanishing of f B £[£] has then to be guaranteed by additional hypotheses, 
possibly in many different ways. Under stronger hypotheses on £ (or on ( , or 
on the boundary conditions which g and g have to satisfy) the set of ADM 
foliations leading to time-conserved quantities with respect to different times 
may be possibly enlarged. Different sets of conditions which guarantee time- 
conservation will be discussed below and we shall compare them to those found 
in the current literature; see also Appendix C for some examples. 

In Section 2 we shall prove that the two action functionals coincide pro- 
vided that the metric g and the background g agree on the boundary dD of 
the region D under consideration. 

In Section 3 the variational principle for York's action functional will be 
reviewed and the definition of quasilocal energy recalled. We shall also briefly 
review the variational principle for the first order covariant Lagrangian. The 
two variational principles are relevant by their own and of course are found to 
be equivalent when g and g are required to agree on the boundary dD , i.e. 
when the action functionals are equivalent. 

In Section 4 a (3 + 1) decomposition of York's action functional will be 
reviewed. The obtained ADM Hamiltonian will be later compared with the 
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Nother conserved quantity of the first order covariant Lagrangian as introduced 
in Section 5 and 6 . 

In Section 5 we shall review the covariant approach to Nother conserved 
quantities and comment on the state of the art in the definition of the energy 
(mass, angular momentum, etc.) in General Relativity. Here most of the ge- 
ometric techniques and bundle framework have been isolated and reviewed in 
Appendix B . 

In Section 6 the Nother theorem will be specialized to define the covari- 
antly conserved quantities of the first order covariant Lagrangian. Such quan- 
tities are involved in the definition of the covariant ADM Hamiltonians (one 
for each spacetime vector field £ ) . They will in fact be compared with the 
ADM Hamiltonian introduced in Section 4 and found to agree (under the 
usual matching conditions required on g and g ) . This result will be obtained 
by comparing the two quantities with respect to the same ADM foliation of the 
region D under consideration. Finally, the relation between covariant ADM 
Hamiltonian and quasilocal energy will also be analysed. The main result is that 
quasilocal energy is covariantly conserved and it is the Nother charge associated 
to the unit normal to the leaves of the ADM foliation. 

In Section 7 we shall discuss different sets of conditions under which con- 
servation in time follows from covariant conservation. 

In Appendix A we collect formulae which are used during the paper to 
translate covariant objects to objects adapted to the ADM foliation and vicev- 
ersa. 

As we already said, Appendix B contains a quick review of the bundle 
framework and Nother theorem at bundle level. 

Finally in Appendix C we present two worked examples which in various 
occasions are quoted throughout the paper. The first one is the computation 
of various Nother charges of the Schwarzschild solution relative to Minkowski 
metric matched on a finite sphere. Conservation in time of various foliations is 
analysed and quasilocal energy is obtained in agreement with previously known 
results (see [3]). The second example is a Kerr-Newman solution matched at 
spatial infinity with the Minkowski metric. The Nother conserved quantities 
are obtained. Such an example is interesting because it does not obey the same 
matching conditions required all over the rest of the paper. It nevertheless 
produces a current which is time-conserved, showing that all the conditions 
discussed along the paper are sufficient but not necessary. 

2. Comparison of the Action Functionals 

From now on we assume dim(M) = 4 , unless explicitly stated. We shall 
here decompose the first order covariant action functional (1.1) along an ADM 
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foliation of D in order to prove that the action functional (1.3) and (1.1) are 
equal if the 4 -metrics g and g are required to coincide on dD . 

Of course, the ADM splitting breaks down the explicit covariance in the 
action allowing a comparison with respect to the same ADM foliation. In partic- 
ular, the boundary term in the covariant action (1.1) splits into a contribution 
on B and a contribution on the lids E io and S tl . 

Let us consider local coordinates (t, r, x A ) adapted both to B and the 
ADM splitting on D . In this coordinate system B has the expression r = 
constant while the leaves T, t are the hypersurfaces of equation t = constant. 
The metric tensor g can be split with respect to the ADM-foliation obtain- 
ing the expression (A. 13) in the Appendix A. Similarly, one can consider the 
foliation of spacetime in the hypersurfaces r = constant and obtain the expres- 
sion (A. 19) in the Appendix A. Analogous expressions can be obtained for the 
reference background metric g . 

Let us evaluate the boundary term of the covariant first order action func- 
tional (1.1) on B and on the lids S fo and E tl , i.e.: 

A B = ~Y K j B >fisr<» dSa =~h ST dSa (2 - 1} 

By using results which are summarized in Appendix A (see equations (A. 18) 
and (A. 23) ) one obtains: 

- Qa f 3 + Y 7 lJ j - {V 1 - V*) (VJ - V*)+ (2.2) 

where V{ is the covariant derivative with respect to the 3 -metric 7^ induced 
on B by the background metric g , while V and V 1 are the radial lapse and 
the radial shift as defined by equation (A. 19) in the Appendix A. We stress that 
no matching condition is required to obtain the above result. The behaviours 
of the metrics g and g are completely unrelated till now. 
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Analogously, on the lids the ADM splitting of the boundary term (2.1) 
gives an extra contribution of the following form: 



Ay H = + — 



1 g^w^u a Vhd 3 x=- -L /" Etl |-2K+ 

+ i? a6 f ^ h ab + Sr h ab ] -^t(N a - N a ) (N b - N b )+ ( 2 - 3 ) 

~ ~M daR (iya " + Iv Da{Na ~ Ra) ) 



where D a is the covariant derivative with respect to the 3 -metric h a b induced 
on E by the background metric g , while N and N a are the lapse and the 
shift of the metric as defined by equation (A. 13) in the Appendix A. Once again 
we stress that no matching condition is required to obtain the result. In general 
(i.e. if no physical requirement about the matching of the dynamical metric and 
the background on the boundary is imposed) the two action functionals (1.1) 
and (1.3) are fairly different. 

However, let us assume that the dynamical metric g and the background g 
coincide on the hypersurface B so that, in particular, they induce the same 3- 
metric on B (i.e. 7^ \ B = 7^ \ B ) and they have the same radial lapse function 
(i.e. V \ B = V \ B ) and radial shift vector (i.e. V 1 \ B = V 1 \ B ). Then, under 
these additional hypotheses, the contribution As reduces to: 

A B = ~J (V7 © - e) d 3 x = ~ J (^0 - v^e) n a ds Q (2.4) 

Analogously, if the metric g and g are required to agree on the lids (i.e. if 
htj = htj , N = N and N l = N l on Yit and ), then the contribution 
(2.3) on the lids reduces to: 

A^=-y (VhK-VfiK^ d 3 x = - J (yg~K - ^ K^j u a ds a (2.5) 

Then the boundary term in the covariant first order action (1.1) can be written 

as: 

Ab+aI'* =-i^( v ^e- v / lG)n Q ds Q +^ H ^K-V§K^u a ds a (2.6) 
and the action functionals (1.1) and (1.3) clearly concide. 
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We stress that this result has been obtained by requiring the aforemen- 
tioned matching conditions between g and g on the complete boundary dD = 
B + S tl — E to of the region D . 

We should also remark that for time-independent solutions the matching 
on the lids cannot be required, because if the two metrics g and g agree 
on a spacelike hypersurface S t , they necessarily agree on the whole region 
D . Also for this reason, the contributions on the lids are never considered in 
applications; this is usually done by restricting to situations in which the lids 
are not present (e.g. by considering a non-compact region D in which to and 
ti are let to tend to — oo and oo , respectively) or are identified. [Identification 
is obtained, e.g., when the solution is time-periodic, as it may happen in the 
Euclidean sector (see [20], [21], [24]) and in approaches based on path-integrals 
for evaluating the grand-canonical partition function or the density of states for 
General Relativity, where the sum over periodic histories has to be considered 
(see [21]). In all those cases the boundary dD is required to have the topology 
<9£ x 5" 1 , i.e. it is assumed to be a single boundary component dD = B .] 

The matching condition of the 4 -metrics g and g is a stronger require- 
ment than the one introduced in [10] . There, only the induced 3 -metrics 7^ 
and 7^ are required to agree on B , where B is let to tend to infinity. Although 
the matching of the 4 -metrics at infinity is not too hard to be implemented 
in applications (see [15], [16], [24]), we are aware that this matching, when 
possible, may become hard to implement in a finite region. Here the matching 
is required to have a direct theoretical comparison between the action function- 
als (1.1) and (1.3) (and between their Hamiltonians and conserved charges). 
In Appendix C we shall explicitly discuss a simple but relevant example of the 
matching of the 4 -metrics in a finite region. 



3. Comparison of the Variational Principles 

We shall here review the variational principles associated to the action 
functionals (1.1) and (1.3) . They are found to agree when the matching 
conditions already discussed in Section 2 are again imposed. We shall also 
recall the definition of quasilocal energy which will be later compared with 
Nother conserved quantities. 

Let us analyse the variation of the action functional (1.3) . It reads as (see 

[7]): 

8I D [g} = ± [ G liv 8g» v ds+ [ (U^5 llJ )d 3 x+ (P ab 5h ab ) d 3 x (3.1) 

JD JB JT,t 
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where G^ = R^ v — | R g^ is the Einstein tensor, while IP- 7 and P ab denote 
the gravitational momenta conjugated to the metrics 7^ and h a b , respectively; 
they are given by: 

n« = -£V7(e 7 y -e«) 

p ab = ±-^fh(Kh ab -K ab ) 

Analogous terms come from the variation 5Id[9~] ■ Notice that, if boundary 
conditions which keep the 3 -metric fixed on the boundary are imposed, namely 
&i\b = and 5/i|s t = 5/i|e 4i = , then the action functional (1.3) is correctly 
extremized by the solutions of Einstein equations = . [We recall that to 
obtain the expression (3.1) the hypothesis of orthogonal boundaries has been 
assumed. Otherwise further terms would appear which have to be cancelled 
through the introduction of a further boundary term into the action (1.3) ; see 
[8], [9], [10].] 

The contribution 5I& = f B (IT J Sjij) d 3 x in (3.1) is considered (see [3]) 
in order to define the quasilocal energy E[g] , i.e. the energy of a region of finite 
spatial extent: 

E,[9] = -/ B ,w d21 (3 ' 3) 

It represents the change of the action in time, where changes in time are gov- 
erned by the lapse N on the boundary B . We stress that the interpretation of 
this quantity as the energy of a gravitating system was suggested in [3] through 
a Hamilton-Jacobi analysis of the action functional. 

Inserting into (3.1) the following decomposition (see [3] for the details): 

*y y = -2^1 5N - 2^pl 5N k + a\a) 5a kh (3.4) 

we obtain: 

E t [g] = 2 I (5^- UiUj) d 2 x=-[ K ^Jad 2 x (3.5) 

JB t \ N J K JB t 

(in the latter equality we have used the definition (3.2) together with the 
expression (A. 11) and (A. 12) in the Appendix A). A contribution analogous 
to (3.5) comes from the background metric g and from the variation Sl^ig] ■ 
Assuming that the induced 3 -metric 7 and 7 are matched on the surface B 
we finally have: 

E t [g,g} = - [ {K-K)^d 2 x (3.6) 

K JB t 

This is the explicit expression of Brown and York quasilocal energy computed 
for the action functional (1.3) . 
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Let us now consider the action functional (1.1) . We can perform the co- 
variant variation with respect to the metrics before choosing any ADM splitting. 
In this way we obtain: 

8A D [g, g] = ± f Sg"" ds - £ / 8g» v ds+ 

J D J D 

r r (3 - 7) 

= h / G^Sg^da-^ / G^8g^ ds+ 

J D J D 

h J ({y/gsT ~ VI <T) Su% v - 8{y/ggn <„) ds a 



+ 



Then, if both g and g are solutions of Einstein equations G^ v = and G^ u = 
, they extremize the action (1.1) , provided that bg^ v = and Su^ v = on 
the boundary dD . 

This is a stronger requirement than the fixing of the metric induced on 
dD . In fact, it amounts to fix the whole 4 -metric restricted on dD (as well 
as the first order derivatives of the reference metric g ) . In any case it is what 
can be done in order to preserve covariance if no additional hypotheses on the 
background are required. 

Nevertheless, we again stress that the action (1.1) is extremized by solu- 
tions of Einstein equations also if we require o~g^ u \ dD = together with the 
matching condition g^lgp = 9^ u \qe> ■ This is obvious from expression (3.7) 
but it also follows from the demonstration carried out in the previous Section: 
under the assumption g^\ dD = 9^v\ dD the action functional (1.1) equals the 
action functional (1.3) . Hence the variation (3.7) turns out to coincide with 
the variation (3.1) . 



4. ADM Hamiltonian 



We shall here review the (3 + 1) Hamiltonian formulation of the action 
functional (1.3) . The ADM Hamiltonian will be later compared with Nother 
charges. 

As it is well known, one of the possible way to give a Hamiltonian de- 
scription of General Relativity is the so-called ADM formulation (see [4], [6], 
[25]). We remark that in this framework the concept of (manifest) general co- 
variance is lost from the beginning owing to the ADM splitting of spacetime 
into space + time. We shall see below how a covariant ADM formulation of the 
problem, deeply related to conserved quantities and Nother theorem, can be 
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formulated. But, first of all, in order to compare known results with ours on 
the subject let us review the standard ADM Hamiltonian formulation of the 
action (1.3) . 

To obtain the Hamiltonian of the system we need to separate out the terms 
in the volume integral ^ J D ^/g~R ds which are pure divergences and then be- 
come boundary integrals. Using the decomposition (A. 13) of the metric we 
have (see [6]): 

^R = ^(n + K^K^ - K 2 ) - 2^~g V„ (Ku» + a") (4.1) 

where TZ is the scalar curvature of the 3 -metric h^ v while a M = u^V ' u u^ is 
the covariant acceleration of the unit normal . Inserting this expression into 
(1.4) we obtain: 

I M=^J d V9(K + K^K^ - K 2 ) ds 

+ \j 1 U ^ Ku>I + a ^^ d3x ~ \f n v( Ku>I + a») v ^d 3 x (4.2) 
+ - / H KVhd 3 x-- f e^d 3 i 



Owing to the properties n M = — 1 and a M = , the integral L H on the 

lids vanishes. Moreover, since we have imposed the condition n fl u ll \ B = of 
orthogonal boundaries, the action functional Io[g] reduces to: 

lD[g} = ^- [ (K + K^K*" -K 2 ) v^ds-- / + (4.3) 

* K JD K JB 

By using formula (A. 12) in the Appendix A we obtain: 

Id[9) = ^[ {n + K^K^-K 2 ) ^gds-- [ IC^d 3 x (4.4) 

ZK JD K JB 

This action may be written in canonical form by inserting the equality: 

K ab K ab -K 2 = ^ \p ah d t h ah - 2P ab D a N b - ^= {2P ab P ab - P 2 ) } (4.5) 
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(see definition (3.2) ) and by removing the term involving the derivatives D a N b 
by an integration by parts. We finally have: 



I D [g] =J dt J^ ( pabd thab -NH- N a H a )d 3 x+ 
~f dt f {~~ IC + Na 2 ~J% Ubaa ^) v ^ d ' 2x 



(4.6) 



where 



n = -^-Vhn+ ^=(2P ab p ab - p 2 ) 

2 « Vh (4.7) 



ab 



n a = -2D b p 

are the Hamiltonian and the momentum constraint, respectively. The Hamilto- 
nian for vacuum General Relativity is thus identified with the term: 

H(g)= [ (NH + N a H a )d 3 x + - [ (NIC - N a K bc n b a ac ) ^ 2 x (4.8) 

(in the latter term we have taken into account that h bc nbO~ ac = ). The Hamil- 
tonian H(g) is the sum of a constrained volume term, which is vanishing when 
computed on a solution, and a boundary term. 

We may repeat the above analysis also for the background action functional 
-Td[<7] obtaining an Hamiltonian H(g) which agrees with (4.8) provided that 
we replace the terms there involved with the corresponding barred ones. The 
total Hamiltonian H(g 7 g) is then given by the difference: 

H{grg) = H{g)-H{g) (4.9) 

When it is evaluated on solutions of field equations, the constrained volume 
terms vanish and (4.9) reduces to the boundary terms: 



H(g, g)~± f (N!C-N a K bc n b a ac ) y^d 

K JB t 

- - I (Nit-N a K bc n b a ac ) Vod : 



(4.10) 



where the symbol ~ denotes equality on-shell. 

Moreover, if the metric g and its relative background g agree on the 
boundary the Hamiltonian simplifies as follows: 

H(g, g) ~ - / {N (/C - K) - N a (K bc - K bc ) n b a ac } ^x (4.11) 
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The mass associated with the time translation = Nu^ + , relative to two 
solutions g and g , is simply defined to be the value of the Hamiltonian (4.10) . 
Clearly the mass H(g,g) of the background is equal to zero. 

Notice that, choosing a Gaussian gauge, i.e. setting N = 1 and N a = 
in (4.11) , we obtain the quasilocal energy (3.6) (see [3]). 

In [26] it was shown that the definition (4.10) agrees with the expressions 
of energy already defined in literature for spacetimes with different asymptotic 
behaviour. 

5. Conserved Quantities 

In this Section we shall analyse conserved quantities associated to the co- 
variant first order Lagrangian (1.1) by the Nother theorem. As we claimed 
above, the Lagrangian (1.1) was originally introduced because it provides a 
simple framework to determine the density of conserved quantities. We shall 
here apply the general framework (see e.g. [14], [15], [16], [27] and references 
quoted therein). 

The Nother theorem is a direct consequence of the general covariance of 
General Relativity (as well as of any other generally covariant field theory). 
It algorithmically defines a Nother conserved current, i.e. a map £[£, cr] which 
associates to any spacetime vector field £ and any field configuration a = (g, g) 
an {n — 1) -form on spacetime M of dimension n which is closed on-shell, 
i.e. when the configuration a is a solution of field equations. For the action 
functional (1.1) we obtain: 



The differential of the (n — 1) -form £[£, a] satisfies the following property: 



where W[f , a] = -(1/2k) (G^ £^ v - G^ v £^ u ) ds is proportional to field 
equations. Consequently, the Nother current £[£,cr] is closed along solutions. 

The equation (5.2) is called weak conservation law for the Nother current 
(5.1) . We stress that equation (5.2) is written on spacetime M , where it 
does not single out a unique Nother current £[£, a] . In fact, if a closed form is 
added to £[£, a] we get another solution of (5.2). This is why Nother's theorem 



- (V%, - SfS?) V a £tgr ~ e R)} ds A 



(5.1) 



(5.2) 
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should be regarded as a claim on (a suitable prolongation of) the configuration 
bundle (see Appendix B), i.e. £ [£, a] has to be regarded as a map £[£] which 
associates a form on (the jet prolongation of) the configuration bundle to any 
spacetime vector field £. As shown in [27] and [28], the Nother current £[£] 
at bundle level is canonically associated to the Lagrangian. Then it can be 
computed along a configuration a to give £[£, a] , i.e. a (n — 1) -form on the 
spacetime M . 

Of course, the Nother currents, as well as the conserved quantities associ- 
ated to them, explicitly depend on boundary terms in the Lagrangian. This is a 
very well known feature in Physics, as it can be simply seen, e.g., in thermody- 
namics. It is in fact well known that boundary terms are related to boundary 
conditions. As we remarked in Section 4 , different boundary terms need dif- 
ferent boundary conditions to keep field equations satisfied by action extremals 
(see [7], [21], [29]). And it is well known that, e.g., for a correct definition 
of energy in thermodynamics, different boundary conditions correspond to dif- 
ferent definitions of energy, such as internal energy, free energy etc. We stress 
that all these energies are true physical energies of thermodinamical systems. 
Which one is to be used in practice is determined by the particular system un- 
der consideration and the boundary conditions we decided to impose. As it is 
physically relevant to notice, we may decide to keep temperature fixed on the 
boundary of a gas box or we may impose adiabatic conditions; this different 
choice corresponds to a different apparatus which selects a different energy flow 
through the boundary so that the boundary conditions are satisfied. We stress 
that this corresponds to an external action on the system which turns out to 
change the physical energy of the system itself. In the covariant first order ap- 
proach to General Relativity something fully analogous holds: the background 
g canonically selects both the boundary conditions and the corresponding en- 
ergy to be used and different choices of the background correspond to different 
physically meaningful definition of energy. 

One can also define (see Appendix B) the superpotential U[£,cr] and the 

reduced current £[£, a] as those currents such that: 

E[£,a] = £[M + M[M (5-3) 

where the reduced current is required to vanish on-shell. 

Once again both £[£, a] and cr] are not uniquely identified by equation 
(5.3) . The superpotential is defined modulo forms which are closed on-shell, 
while the reduced current is defined modulo forms vanishing along solutions. 
The decomposition (5.3) of Nother's current is again well-defined only at bun- 
dle level, where a decomposition algorithm can be constructed (see [28], [27]). 
The bundle superpotential W[£] and the reduced current £[£] are canonically 
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and globally denned and then they are computed along the configuration a to 
give U[^,a] and £[£, cr] , respectively. 

The conserved quantity in a region O is defined as: 

Qn[Z,*]= [ £[M= [ [ U[&<r]= [ U[£,v] (5.4) 

JO. JQ JdQ JdQ 

where, in the last equality, a is assumed to be a solution so that £[£, cr] = . 

In the case of the first order covariant action functional (1.1) we obtain in 
particular: 



= -\y/g 9^ C-Vti r x £ 



K 



ds a/3 

(5.5) 



where g^ v (as well as Christoffel's symbols ) has to be regarded as local 
coordinates on (the jet prolongation of) the configuration bundle Lor(M) of all 
Lorentzian metrics on M . They become the metric components, i.e. functions 
of spacetime point x G M , only when they are calculated along a configuration 
a(x) = (g(x),g(x)) . 

The conserved quantity <3n[£, cr] depends on the dynamical metric g and 
on the reference background metric g . It has to be interpreted as the relative 
conserved quantity of g with respect to g . For example, if the energy is con- 
sidered (by choosing the vector field £ in a suitable way, see below), Qnf^c] 
represents the amount of energy which is necessary to pass from the spacetime 
(M, g) to the spacetime (M, g) . As it is physically reasonable, this energy 
can be infinite in principle. Of course, if the metrics are "near" the energy 
between them can be expected to be finite, though we cannot even be sure 
that we can "continuously" join two very different metrics. In general we can 
expect conserved quantities to classify accessibility classes of metrics, though 
we do not want here to enter the problem of providing the set of sections of the 
configuration bundle with a topology or a differentiable structure. We certainly 
expect that if a = (g, g) , i.e. if we set g = g , the conserved quantity Qn[£, cr] 
vanishes. This condition is satisfied by prescription (5.4) because of the form 
of the superpotential (5.5) . 

We remark that the conserved quantities defined by (5.4) are covariantly 
conserved, meaning that their flows through boundaries dD of 4 -dimensional 
regions D in spacetime M identically vanish. If one chooses g to be a solution 
asymptotically flat according to one of the current definitions (e.g. the Kerr- 
Newman solution), g to be the flat reference background (which matches the 
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dynamical metric at infinity where "infinity" is prescribed by the definition of 
asymptotic flatness) and O to be a spacelike hypersurface in M , then cr] 
reproduces the expected value for mass (by chosing £ = dt , i.e. the vector field 
which corresponds to asymptotic time translation) and angular momentum (by 
choosing £ = — , i.e. the vector field which corresponds to asymptotic rota- 
tion; see e.g. [11], [15], [27]). Furthermore, the same results are achieved for 
non-asymptotically flat solutions by choosing suitable reference backgrounds. 
Examples are the (2 + 1) BTZ solution (which is asymptotically anti-de-Sitter, 
see [16]), the Euclidean Taub-Bolt solution (which is asymptotically locally 
flat, see [30]). In addition, the same techniques are used successfully in gauge- 
natural theories, i.e. when the field theory owns both covariance and gauge 
invariance (e.g. BCEA theory, Einstein-Maxwell theory; see [27], [31]). 

We stress that the conserved quantities Qn[£> cr] associated to the covari- 
ant first order action principle (1.1) are not affected by the anomalous factor 
problems as the ones associated to the standard Hilbert-Einstein Lagrangian 
(see [14]). As is well known, in fact, the Komar superpotential (see [32]): 

U K = ^V^ a ds aP (5.6) 

which is the superpotential associated to the Hilbert-Einstein Lagrangian, when 
computed, for example, on the Kerr-Newman solution, produces the correct 
angular momentun, but just one-half of the expected mass. Of course one could 
postulate the mass to be associated to the vector field £' = 2dt ■ However, 
the factor appears to depend on the particular solution under investigation, 
thus such a prescription seems to be incorrect as well as unmotivated. This is 
usually interpreted as a hint of a correction needed in the definition of conserved 
quantities. The correction can be obtained by the ADM techniques (see [4]) by 
restricting to asymptotically flat solutions and by choosing an ADM foliation of 
spacetime M . As shown in [33], [34], [35] and Section 6 , the ADM formalism 
appears as a particular case of the technique exposed above, which furthermore 
applies to much more general situations. 

6. Covariant ADM Hamiltonian 

We shall here specialize the general framework introduced in Section 5 to 
the first order action functional (1.3) . We thence obtain a map Q[£] which 
associate to each spacetime vector field a covariantly conserved quantity, called 
the covariant ADM Hamiltonian. It will be compared with the standard ADM 
Hamiltonian introduced in Section 4 and with the quasilocal energy defined in 
Section 3 . 
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Let us consider the covariant conserved quantity QeJ^ct] for the first 
order action functional (1.1) in the domain T, t and relative to a vector field 
£ and a section a = (g, g) . As we have already outlined in the previous 
Section, the quantity Qs t [£, cr] is defined to be the integral of the superpotential 
(5.5) on the 2 -dimensional surfaces B t (see Fig. 1 and Appendix A for the 
notation), i.e. 



In order to simplify the ADM decomposition of the expression (6.1) so to be 
able to compare the results obtained with the standard ones of Section 4 and 
[3] let us assume, as usual, that the metrics g and g are matched on the 
hypersurface B and that the boundaries are orthogonal (i.e. u^n^g = ). 

[We stress that under our viewpoint the matching condition between g 
and g is unessential, since Nother currents are covariantly conserved. One may 
consider the second example analysed in Appendix C where the Kerr solution 
is studied and its mass inside the finite R-sphere is obtained with respect to a 
flat background matched at infinity. In this Section we require the matching on 
B in order to compare the Nother charges expression with the aforementioned 
standard (3 + 1) Hamiltonian and quasilocal energy] 

Let us also assume that the vector field £ is tangent to the hypersurface 
B (i.e. e n v\i3 ~ )• First of all let us consider the first contribution Qs t [£, g] 
into (6.1) , i.e. the integral of the Komar superpotential. It may be rewritten 
as: 

Q*A^9] = ^f (upn a -u a np)g^V^ a ^d 2 x (6.3) 

*K J B t 

where y/acPx, is the volume element on B t . We can manipulate algebraically 
the latter expression in the following way: 



QTM, cr] = Q-£ t [£, g] + Q Et [£, g, g] + Q St [£, g] 



(6.1) 



where: 





(6.2) 






(6.4) 



2k 



1 



£« V M n Q - n a £^g^ a } ^d 2 x 



JB t 
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By means of the identity n a £^g fia = £^n^ — g^ a £^n a we obtain: 

Taking formula (A. 9) repeatedly into account together with the condition of 
orthogonal boundaries u^n^\ B = and the condition £ M n M | e = we finally 
obtain 

Qv t [Z,g] = - I {e^^T} V^d 2 x+ ^ / u a £sn a ^d 2 x (6.5) 

« JB t 2k JB t 

A similar expression may be found for the third contribution Qj; t [£,g] into 
formula (6.2) , i.e. the Komar contribution of the matched bakground g : 

Qv t [S,9] = — [ {e Ma ^£ Q } V^d 2 x-±- f u a £^n a ^d 2 x (6.6) 

K J B t ZK JB t 

It now remains to calculate the second contribution Qs t [£, fiS i n t° formula 
(6.2) . We stress that this is the contribution arising from the boundary term 
into the action functional (1.1) . It can be written as: 



<5s t [£, g,g] = J {upn a - u a np)£, a w& v g» u ^fad 2 x 



(6.7) 



We remind that in the radial ADM decomposition (A. 19) of the metric we 
have np dx 13 = V dr where V is the radial lapse. Hence, by making use of the 
expression (A. 23) in the Appendix A we obtain: 

Qv t [Z,9,9\ = -- I u a C(0-e)V^d 2 x (6.8) 

K JB t 

[Notice that only the projection u a £ a of the vector field £ Q along the timelike 
normal u a gives a contribution to the term Q St [£, g, g] . This is the reason 
why the boundary term into the action (1.1) allows to correct the anomalous 
factor of the Komar superpotential which, as we said above, appears in the 
computation of mass while it does not enter in the computation of angular 
momentum.! 
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The conserved quantity in the region T, t relative to the infinitesimal gen- 
erator of spacetime symmetries £ is given by the sum of (6.5) , (6.6) and 
(6.8) : 



K JB t 
% K JB t 



(6.9) 



This latter formula may be recasted in a form which is better suited to be 
analysed. Because of IC^u^ = and a^u^ = , from formula (A. 11) in the 
Appendix A it follows in fact: 

© M a r = ~u a T a» - < r K pP (6.10) 
Inserting this latter expression together with (A. 12) into (6.9) we obtain: 



+ — / u a (£tn a - £tn a )Vo z d 2 x 
2« 7s t 



(6.11) 



Let us stress that, until now, no assumption has been made on the vector field 
£ , apart from the requirement £ M ?v|g = . 

An easy computation shows that the difference u a (£^n a — £^ n a ) is al- 
ways zero if the metrics g and g are matched on B . 

On the contrary the terms u a £^ n a and u a £% n a appearing in the last 
contribution to (6.11) separately disappear if we choose £ to be tangent to the 
2 -surfaces Bt , that is ^ fJ, u fx \ B = , or also if we choose £ to be the time-like 
vector field dt , i.e. £ a = Nu a + N a . In both cases the flow of the vector field 
£ maps each hypersurface E t into itself or, respectively, into another surface 
S t / . Since the vector field n a is tangent to each E t it turns out that also 
£^n a is tangent to Et and then u a £^n OL = in these cases. 

Hence if we specialize formula (6.11) for the vector field £ = dt we obtain 
the covariant conserved quantity which we call the Hamiltonian of the system. 
It is given by the expression: 

Q™[d u <j] = - [ {N{K-K)-N a {K bc -K bc )n h Gac} V^d 2 x (6.12) 

K JB t 
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and it coincides exactly with the expression of the (3 + 1) Hamiltonian (4.11) . 

Let us notice that this definition of Hamiltonian can be correctly considered 
as the definition of a covariant ADM formulation (see [33], [34], [35]). In fact 
it does not require, a priori, a ( 3 + 1 ) decomposition of spacetime. We stress 
that in the covariant ADM approach, the Hamiltonian, or energy, contained in 
a 3 -dimensional region O and relative to a solution a , is defined by (5.4) as 
a Nother conserved quantity: 

Q3*[£,<7]= / U%d\ (6.13) 

JdQ 

This is a well-posed definition of Hamiltonian provided only that the non- 
vanishing vector field £ be transverse to the hypersurface O . Hence, by consid- 
ering the parameter of the flow of £ as the "time" parameter and transporting 
Q along the flow of £ we obtain a world tube foliated by hypersurfaces diffeo- 
morphic to O . In this covariant context, rather then starting from a preferred 
local foliation into hypersurfaces, the starting point is a non-vanishing vector 
field the flow of which defines the local time, i.e. the flow of evolution. Then, 
by specializing the definition (6.13) to the ADM foliation depicted in Fig. 1, 
under the additional assumptions of orthogonal boundaries and of the matching 
between the metric and its background, the covariant Hamiltonian Qj?*[dt, o~] 
exactly coincides with the standard Hamiltonian H{g, g) derived from a ( 3 + 1 ) 
splitting of the York action functional (see equation (4.11) ). 

Another relevant Nother conserved quantity is obtained by specializing 
formula (6.11) to a vector field £ a = N a tangent to the 2 -surfaces Bt , i.e. 
N a u ol \ Bt = and N a = N p . Because of the vanishing of the first and third 
term in the right hand side of (6.11) we obtain: 

Q%[N a , *] = -- [ N a nP(K af) - K af3 ) ^d 2 x (6.14) 

In asymptotically flat spacetimes when £ corresponds to a rotation at spatial 
infinity, the Nother charge (6.14) may be taken as the definition af angular 
momentum. 

The last Nother charge we consider is the one relative to the unit vector 
field £ = u normal to the leaves of the ADM foliation. From (6.11) we obtain: 

Q£*[u, a] = - f{K-K) V^cfx (6.15) 

K JB t 

We observe that it agrees with the definition (3.6) of quasilocal energy, i.e. it 
is the value of the Hamiltonian (6.12) with N\ Bt = 1 and N a \ Bt = . In 
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the aforementioned hypotheses, quasilocal energy may then be considered as a 
Nother charge. 

7. Time Conservation 

We shall here discuss two different sets of sufficient conditions for time- 
conservation. The quantities (6.12) , (6.14) and (6.15) , are all covariantly 
conserved quantities independently on the hypothesis that £ is a Killing vec- 
tor field or not. In fact they have been defined by means of Nother theorem 
through a construction which relies only on the covariant nature of the La- 
grangian. Hence, on a solution of field equations, the covariant conservation 
law d^S^l^a] = always holds for the Nother current £[£,cr] and for all 
vector fields £ . This property, together with the property of existence of su- 
perpotentials for any natural theory, has allowed us to define the covariantly 
conserved Nother charges Qn ot [£, cr] (see (5.4) ). On the contrary the charges 
<5^ ot [C,c] are conserved in "time" if they do not depend on the chosen hyper- 
surface , i.e. if Qj?*[£, cr] = Q™ [£, cr] . This is a stronger condition that has 

to be supported by additional requirements. If B is the 3 -dimensional region 
such that E t / — E t + B is the boundary of a region D , from the conservation 
law <i M £ M [£, a] = we obtain a time-conserved quantity if jgE^ds^ = , i.e. 
if the net flow of the Nother current S through the hypersurface B vanishes. 
A stronger condition amounts to require the integrand to be equal to zero on 
B , i.e. £ M ?v|g = . In this case Q Tot [£, cr] is conserved not only with respect 
to the given foliation in hypersurfaces E t but it is time-conserved with respect 
to the time of any foliation of the region D . 

For the action functional (1.1) the Nother current (5.1) may be rewritten 

as: 

£ a [M = ^ {(Vgg^ - vW)^„ - ^(VsO - £ Q £} (7.1) 

The corresponding quantity £^n^\ B will be equal to zero if some condition 
is verified. We do not explicitly known a set of necessary requirements for 
the occurence of this situation; we can only provide two examples of sufficient 
conditions. 

We may require, as a first example, that the following properties hold true: 

A) the vector field £ is a Killing vector field for the metric, i.e. £^g^ v = ; 

B) the vector field £ is a symmetry for the background in the sense that 

C) £ is tangent to the boundary B , i.e. = . 
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These three requirements together ensure that S^n^\ B = . These prop- 
erties are clearly satisfied if we are dealing with the Killing vector fields of an 
asymptotically flat stationary solution and we choose the flat metric as a back- 
ground. Nevertheless, we stress that the time-conserved quantities Qj?*[£, cr] = 
J B U[£, cr] can be computed on a finite region, i.e. it is not necessary that B t 
is identified with spatial infinity. Moreover, we do not have explicitly required 
here the matching between the metrics on B (see the example of the Kerr 
metric in Appendix C). 

Another set of sufficient conditions may be imposed to fulfill the condition 
E^n^\ B = 0. They closely resemble the ones of [3]. We may require the 
matching of the metrics on the boundary B (so that the first term in (7.1) 
vanishes) and we again require £ to be tangent to the boundary B : £^ = 
(in order to make the third term in (7.1) vanishing when contracted with the 
normal n a ) . We are left with the term: 

J S a da a = ~ J My/99 1 ") ds « 

Owing to the matching requirement g^ u \^ = 9hAb ^ ne l & tter expression may 
be recasted into the equivalent form: 

/ £ a ds a = [ (U ij -U ij )£ aij d 3 x 

JB JB 

It vanishes if we require £ to be a Killing vector of the boundary 3 -metric: 
£^ij = T>i^j + T>i^j = . Hence, also in this latter situation, we obtain 
time-conserved quantities Q^[^,cr] , for the time t of any foliation of D in 
hypersurfaces E t . 

8. Conclusion and Perspectives 

We proved that once suitable matching conditions are required (i.e. the 
4 -metrics g and g are required to agree on the boundary of the region under 
consideration) the two action functionals (1.1) and (1.3) agree. Consequently 
the action functional (1.3) may be considered as the ADM counterpart of the 
covariant action functional (1.1) . 

The second important result achieved here is the characterization of the 
quasilocal energy as the Nother charge associated to the (timelike) unit vector 
normal to the leaves of the ADM foliation. 

These seem to be new results which should enable us to extend the analysis 
further ahead to the prescription for the entropy in General Relativity. In fact, 
the quasilocal energy as well as the action functional (1.3) appeared also as 
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the starting point of a statistically-oriented approach to black hole entropy (see 
[21], [22], [24], [26], [36], [37] and references quoted therein). A different 
approach to black hole entropy based on Nother approach (see [16], [27], [30], 
[38] and references quoted therein) may be found in literature. In view of 
the present comparison between the covariant first order approach (which the 
Nother approach is based on) and the York's action functional (which quasilocal 
energy is based on) as well as between the Nother charges and the quasilocal 
energy themselves, we believe that the two different approaches to entropy can 
be now successfully compared (see also [39], [40]). 

Another interesting perpective is to extend the present comparison to the 
more general situation of non-orthogonal boundaries which sometimes appeared 
in the literature (see [8] , [9] , [10] ) . It would be of some interest to know whether 
non-orthogonality of the boundaries in the ADM decomposition of the covariant 
action functional (1.1) exactly produces the additional boundary terms which 
are derived for the (modified) York's action (see [9]), as we guess to be true. 
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9. Appendix A: Notation 

In order to make this paper self-contained we here briefly summarize some 
of the formulae and expressions which are used throughout the paper. We follow 
the convention and notation adopted in [3] and we also refer the reader to Fig. 
1 for notation. 

We assume the hypersurfaces E$ to be spacelike and we assume the hy- 
persurface B to be timelike. The metric induced on the hypersurfaces E t 
may be written as: 

hpv = g^u + u^u u (A.l) 
while for the metric 7^„ on the hypersurface B we have: 

liiv = 9iiv - % n v (A.2) 

(in the sequel Greek indices are always raised and lowered with the 4 - dimen- 
sional metric) . The two vectors u and n denote the future directed unit normal 
to Ef and the outward pointing unit normal to the hypersurfaces B , respec- 
tively. They satisfy the normalization relations u^u^ = — 1 and n^n M = 1 , 
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respectively. Any spacetime tensor may be projected onto the hypersurfaces E t 
by means of the projection tensor: 

ht = 6£ + u^u v (A3) 

Any tensorial object may be also projected onto B with the projection tensor 

tf = 8$-n»n v (A4) 

The 2 -metric on the boundaries B t is given by: 

cv = + u fl u u = h^ u - n M n u ( A5) 

and the respective projection tensor is a£ = g^ p a pv . 

The extrinsic curvatures K^ v of E ( in M, Mi , of B in M and JC pu 
of i?t embedded in T, t are defined, respectively, as follows: 

©^ = -7^V a n, (A6) 

where _D a denotes the covariant derivative on S t compatible with the metric 
h. 

The extrinsic curvature .fT^j, is a symmetric tensor on S t , i.e. it satisfies 
the conditions K pv u^ = , K pv h^ = K pv . Instead the extrinsic curvature 
is a symmetric tensor on B : O^n^ = , Q^lp = ® P v > while the extrinsic 
curvature /C M „ is a symmetric tensor on B t , i.e. JC pu u^ = JC pu n p = . 

We also denote by: 

a" = u»V u u v 

(A7) 

b u = n^V p n v 

the (covariant) accelerations of the two normals w M and , respectively. They 
satisfy the orthogonality properties: a M = and n M 6^ = . 
By making use of the property (A4) , we obtain: 

V,n^ = V a n^ = 7« V a n^ + n a n„V Q ?^ (A8) 

Taking into account definitions (A6) and (A7) we have: 

V un v = _Qv + n u b» (A9) 
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Performing calculations in the same manner for the vector , it is easy to 
check the analogous relation: 

V v vP = -K£ -u v aP (A10) 

Moreover, projecting the indices of the extrinsic curvature Mi , normally and 
tangentially to the hypersurfaces E t we obtain the useful formula (see [3] for 
detailed computations) : 

Qfiu = k-iiu + u v n a a a + 2a^u v) n p K af3 (All) 

[We remind that this relation is true only under the assumption of orthogonal 
boundaries, i.e. n^\ B = .] Contracting the latter expression with g^ v we 
also easily obtain: 

e = /C-n a a Q (A12) 

Let us now consider the ADM decomposition of the metric: 

g = g^dx^ (g) dx v =(h fliy — u^u^dx^ ® dx v = 

( A 13) 

= — N 2 dt 2 + h a b(dx a + N a dt) ® (dx b + N b dt) 

The coordinate system (x M ) = (£, x a ) = (i, r, x A ) is adapted both to B and to 
the ADM foliation on D . The surfaces E$ are surfaces of constant t while B 
is the hypersurface of constant r . In other words, indices a, 6, c, • • • run from 
1 to 3 and denote indices on the spacelike hypersurfaces E t , while indices 
A, B, C, • • • run from 2 to 3 and they instead denote indices on the boundary 
B t of E t . Hence, tensors on E t are labelled by early Roman letters a, 6, ... . 
When they are considered as tensors on spacetime M the same tensors are 
instead denoted by Greek letters. For example, the extrinsic curvature K in 
(A6) can be denoted as K^ v or K ab , according to notational convenience. 
The unit normal is given by (it M ) = (—AT, 0,0,0) while the timelike co- 

ordinate vector field d = d/dt reads as d = Nu + N , being N the lapse 
function and N = (A/>) = (0, N a ) the spatial shift vector: N ■ u = = 

(we remind that if we also assume the vector field do be tangent to the 3- 
dimensional hypersurface B , the orthogonal boundaries condition reads as 

N ■ n = n^lg = . Nevertheless this latter hypothesis will be not relevant 

for the computations in the rest of this Appendix) . The extrinsic curvature 
and the acceleration defined in (A6) and (A7) , respectively, read as: 

K ab = -NF° ab = [-d h ab + D a N b + D b N a ] 

\ {AAA) 
r) b N 

(O = (0y) = (0,— ) 
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(Roman indices are here raised and lowered with the 3 -metric h ab ). 
The Levi-Civita connection coefficients are given by: 



Oa 



N 



(d a N-N b K ba ) 



r° 



ab 



00 



-pa 
1 bO 



ra 
be 



_Kab 

N 

d Q N N b N a N b 

— + JT dbN -— Kab 

N a N c N a 
-NK% + K cb + D b N a - —d b N 



3pa 



N a 

be + ~Jf K bc 



(A.15) 



(where 3 T^ c denotes the Levi-Civita connection of the 3 -metric h ab ). A sim- 
ilar splitting may be performed with the background metric g ; we obviously 
obtain the same relations by replacing the objects involved with the corre- 
sponding barred ones. By means of formulae (A15) one may easily compute 
the quantities u^p = — S^ a T e ^ e . For example, one can verify the following: 



(A16) 



and 



' (N a -N a )(N b -N b )- ^ (^- T h ab + 



N 2 N 
1 



ab 



N \N 

N b 



+ 



I N b N b 

v y2 - N a )d a N + — 2 d a N« + ^ 3 i- - — 2 *r: b 



which together give the expression 



(A17) 



(A18) 



The latter term is involved in the contribution on the lids of the boundary 
term in the action functional (1.1) . To evaluate instead the action functional 
contribution on the hypersurface B we have to make use of the adapted splitting 
of the metric: 



9 = 9^dx^ <g> dx u =(7^ + n^n u )dx^ <g> dx u = 

=V 2 dr 2 + ^ ij (dx i + V 1 dr) <g> (dx J + V j dr) 



(A.19) 
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where (x l ) = (t,x A ) . Middle Roman letters i,j,k,--- denote indices on the 
timelike hypersurface B while x A (A = 2,3) denote again the coordinates 
over the two dimensional surfaces B t . The function V is the radial lapse 
while V 1 is the radial shift. Hence the unit, outward pointing, radial normal n 

reads as n = (1/V)[d r — V l di] . The extrinsic curvature 0^ of the "cylinder" 
B = {r constant} is given by: 



e 



VY r - 



2V 



[-drlij + ViVj + VjVi] 



(A20) 



where T>i denotes the covariant derivative on B induced by the Levi-Civita 
connection 3r *fc of the 3 -metric 7^ . The coefficients of the 4 -dimensional 
Levi-Civita connection can now be decomposed as: 



r r . = 

rt 


±_ {div+v j Qji) 




r r . = 


V 




T r = 

rr 


1 V* 

-d r V + —d t V + 


v i v j 
v e « 


r*. = 


v i v h 
-ve) - — e kj 


+ VjV* - 


r* — 
1 jk - 


v i 

3 F* ft 

ijk - y jk 





{A.21) 



V 



The latter expressions are not simply obtained from (A. 15) by exchanging 
tensors on £ t with the corresponding tensors on B . Because the metric h a \, 
and 7y have different signatures, a change of sign may appear in some terms 
of (A.21) if compared with the decomposition (A. 15) . 

By means of (A.21) , the following expression are then easily computed: 







^x, = 


V 1 




% 




1 


+ 


vv 2 



1 

V 2 



1 (V 1 -V l )(V J -V J )+ 1 



V 2 V 



V 



V \v 



_.. V 
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+ (A.22) 



1 _ yj yj 
ty 1 _ y l )f) V —8 V l - — 3 H ■ 4- — 3 H • 



Expressions (A.22) give the expression 

g^w^n a = g^(u% - t$> a = Vg^(u^ - u^) 
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(A23) 



which is the contribution on B of the boundary term in the action functional 
(1.1). 

We stress that in (A. 18) and (A. 23) computations are performed without 
any hypothesis of orthogonal boundaries and without requiring any matching 
conditions between the metric g and its background g . 

10. Appendix B: Bundle Formalism for Variational Calculus 

We hereafter briefly recall the bundle framework for variational calculus 
and summarize how one can algorithmically construct the Nother conserved 
currents and the superpotentials out of the Lagrangian. 

Let C = (C, M, 7r; F) be a bundle and (x^;y l ) be fibered coordinates 
(relative to a trivialization chosen on C ). In fibered coordinates the projection 
reads as n : (x M ;y l ) i— > x M . A local section defined on U C M is a map 
a : U — > C such that no a = lu , i.e. it is locally given by a : x^ i— > (a^; y l (x)) . 
A local section is thence associated to a map a 1 : i— > y l (x) . When U = M 
the section is called a global section on C. In variational calculus y l (x) are 
identified with the values of dynamical fields at the point x G M . The bundle 
C is called the configuration bundle. 

Starting from the bundle n : C — > M we can define the bundle n k : 
J k C — > M which is called the -order jet prolongation of C . A point in J k C 
is denoted by j k a and it is an equivalence class of local sections (defined in a 
neighbourhood of x G M ) having contact of order fe at x , i.e. having the same 
Taylor expansion at x up to order . Points in J k C are then parametrized 
by the value of a section a and its partial derivatives at the point x G M up 
to order /c included. Consequently, if (x^;y l ) are fibered "coordinates" on C, 
(x^ 1 ; y l , y l , . . . , ) are fibered coordinates on J k C . We remark that the 

"coordinates" Mh (h < k) are symmetric in the lower indices. If a is a 

section of C , we can prolong it to a section _7 fc cr of J k C which is defined by 
j k a(x) = j k a . If a : x^ i— > (x M ;y l (x)) , the /c -order prolongation j fc cr is given 
by 

/ a . x m _ ^ y i (a;)j . . . ? ^...^^(x)) (B.l) 

If a section cr is identified with a field configuration, j k a describes fields and 
their derivatives up to order k . 

Let us denote by $ : C — > C" a bundle morphism projecting onto a 
diffeomorphism : M — > M' . By definition of bundle morphism it preserves 
the fibers (i.e. 7r ; o $ = o 7r ) and it is locally given by 

f x' = <t>{x) 
ly =$(x,y) 
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We can prolong it to a bundle morphism j k & : J k C — > J k C defined by the 
following (composition preserving) rule 

j'^J^-jJw^^f 1 ] (5-3) 

(see [41], [42] for further details). Analogously, for any infinitesimal generator 
of automorphisms, i.e. any projectable vector field S = ^(x) <9 M + y) <9; on 
C , we can define the prolongation j k E = £ M <9 M +<f di+Qd? + . . •+^ 1 ...^ fc 9f 1 -' ifc 
(here we set = , = ^ , . . . , d^"^ k = ^7 — — for the natural basis 

induced in the tangent space TJ k C by fibered coordinates). It is recursively 
given by 

4 = d M e - = d M (e - y m + y ^ 



P = d P —v 1 8 f u = d f£* - iff 1 ') + ?/ £^ 

(B.4) 

where d M denotes the fotaZ derivative at any order and it is (locally) defined by 

d M = / , + yj l < + y£ M $' + ... (B.5) 

We define a variational principle of order k to be a pair (C, L) where C = 
(C, M, 7r; F) is a bundle called the configuration bundle and L : J k C — > A n (M) 
is a (vertical) bundle morphism called the Lagrangian ( of order k ). Here and 
below Ah{M) (h < n) is the bundle of h -forms on M and n = dim(M) . 
Locally a Lagrangian reads as L = £(x^, y l , y 1 ^, . . . , y^ 1 ... Atfc ) ds (where ds = 
da; 1 A ... A dx n ); £ is the Lagrangian density and it depends on the spacetime 
point x^ , on fields y l and on derivatives of fields (j/^, . . . , y^ tl ... Mfe ) up to a 
finite order . The Lagrangian L associates to any configuration a an n - 
form L o j' fc cr = (Co j k a) ds over the spacetime M . 

In General Relativity, the configuration bundle is the bundle C = Lor(M) 
of Lorentzian metrics on the spacetime M . It has local fibered coordinates 
(x^^g^v) . The Hilbert Lagrangian L = (l/2«) yfgR ds is second order, thus 
it is interpreted as a map L : J 2 Lor(M) — > A n (M) . The second order jet 
bundle J 2 Lor(M) has natural fibered coordinates (x M , g^, g^x, g^Xa) or 7 
equivalently, it can be parametrized by the more convenient set of variables 
(a; M ,^ Ml/ ,r^,d CT r^) where denote the Christoffel symbols of the metric 
9- 



32 



We stress that in this framework, the Lagrangian is a map between finite- 
dimensional manifolds. It induces the action functional: 



Ad(*) = [ 



Loj k a 



(B.6) 



J D 



where D C M is a region in spacetime M , i.e. a compact submanifold (of 
dimension n ) with a boundary dD which is again a compact submanifold (of 
dimension n — 1 ). 

Let now C = (C, M, tt; F) be a fiber bundle. We can define the sub- 
bundle r v : V(C) — > C of the tangent bundle r : TC — > C , which turns out 
to be a vector sub-bundle of TC . A vector v G T p C belongs to V P C if and 
only if T p tt(v) = ; in this case we say that v is a vertical vector at p e C . 
Accordingly, V(C) is defined to be the kernel ker(T7r) of the fiberwise linear 
map Tn : TC — > TM . A vertical vector is "tangent to the fiber" F and its 
local expression (using fibered coordinates {x^;y % ) on C ) is v = v % di where 
<9j = is a local pointwise basis of vertical vectors. A section of the bundle 
V(C) — > C is called a vertical vector field. When C is the configuration bundle 
of a field theory, a vertical vector field X = (5y l ) di represents a deformation 
of the dynamical fields. We can also consider the dual vector bundle V*(C) of 
pointwise linear functionals on vertical vectors, which is called the dual vertical 
bundle. A pointwise basis dy l of V*(C) can be defined to be the dual basis 
to di , i.e. dy'idj) = 8) . Notice that, even if V{B) C T(C) , V*(C) is not a 
sub-bundle of T*(C) . 

We stress that the dual vertical bundle can be defined for any bundle; we 
can consider the dual vertical bundle V*(C) of the configuration bundle as well 
as the vector bundle V*(J k C) of its prolongations. The importance of these 
objects is that they allow to translate the variation of the action functional 
into the variation of the Lagrangian morphism. For any Lagrangian L : J k C — > 
A n (M) we can in fact define a bundle morphism 5L : J k C -> V* (J k C)®A n (M) 
such that for all vertical vector fields X : C — > V(C) we have: 



where $ s : C — > C is the (vertical) flow of X and < > denotes the 
canonical duality between V*(J k C) and V(J k C) given by evaluation. Locally, 
the morphism SL is given by 




s=0 




8L= [pidy'+p? dyi + .-.+ph 



dy^...^] ®ds 



(B.8) 
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where we defined the conjugate momenta (pi = diC, pf = d^C, . . . , p^---^ k = 
df 1 ' Mfe £) to be the derivatives of the Lagrangian density C with respect to 
the "coordinates" (y\ y 1 ^ . . . , ^...mJ > respectively. 

The variation of the action along a deformation X e V(C) is thence given 

by 

<S x A D (a) = / < <5L | j fc X > oj k a (5.9) 

Field equations now follow by Hamilton principle, i.e. requiring the action to 
be stationary along any deformation around classical solutions. In other words, 
we say that a configuration a : M — > C is a classical solution if, for all regions 
D and for all deformations X such that j k ~ 1 X\gD = (i.e. such that the 
deformation X vanishes together with its derivatives up to order k — 1 on the 
boundary dD of the region D ) we have: 



d x A D (a)=0 (B.10) 

Specializing to the Lagrangian (1.1) , which is first order in g and second order 
in g , the deformation X = X g + Xg = (Sg fliy )d fll/ + (Sg^d^ splits into a 
deformation X g of the dynamical metric which is required to vanish on 
the boundary dD (i.e. X g \gD = 0) and a deformation X g of the background 
metric g^ u which is required to vanish together with its first order derivatives 
(i.e. j 1 X g \dD = ) on the boundary dD . 

One can prove that there exist a (unique) global bundle morphism E : 
J 2k C — > V*(C) <S> A n (M) (called the Euler-Lagrange morphism) and (a family 
of) global bundle morphisms F : J 2k - x C -> V*(J fc_1 C) ® A n _ x (M) (called 
the Poincare-Cartan morphisms) such that the so-called first variation formula 
holds true, i.e.: 



< 5L | fX > oj k a =< F | X > oj 2k a + d(<F\ f~ x X > oj 2fc "V) (£.11) 

Since in the Hamilton principle we consider deformations such that j k ~ 1 X\go = 
the Poincare-Cartan contribution vanishes when integrated in (B.9) and we 
get field equations F o j 2k a = , i.e. Euler-Lagrange field equations. 
If, for example, we consider k = 2 we obtain 

E=\pi- cWf + ^vPTW ® ds 
F=M-^)W+p^Ayl]®ds il 

A bundle C = (C, M, 7r; F) is called natural if there exists a natural lift 
to C of spacetime diffeomorphisms, i.e. if one can associate to any spacetime 
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diffeomorphism : M — > M a (unique) bundle automorphism $ : C — > C 
so to preserve compositions. As a consequence, on natural bundles one can 
canonically and naturally lift spacetime vector fields £ = ^(x) 9^ to bundle 

vector fields £ = £ > ^(x)d fl +£ > l (x,y)di . A field theory is called a natural theory if 
its configuration bundle C is a natural bundle and, in addition, any spacetime 
diffeomorphism is a Lagrangian symmetry, i.e. we have 

< 5L | j k £ i a > oj k a = £^(L o j k a) (5.13) 

where £^(L o j k a) = (d o + o d)(L o j u) is the Lze derivative of the 

Lagrangian and i^cr = TV o £ — £ o a is the Lze derivative of the section a 
along a spacetime vector field £ . Locally we have 

£<■* = (£" d^ l (x) - f(x, y(x))) (5.14) 

Notice that i^cr is a vertical vector field on the section a so that the l.h.s. of 
formula (5.13) is meaningful. We also recall that equation (5.13) is equivalent 
to the requirement that the so-called Poincare-Cartan form is invariant under 
the flow of £ (see [2]). 

The global Nother theorem is now easily obtained. Let L be a natural 
Lagrangian; inserting (5.11) into eqation (5.13) we obtain 

< F | £^a > + d < F | j k ~ 1 £^a >= d (i c L) 

i.e. 

d(< F | j k ~ 1 £^a > -i € L) = - < ^ | ^£7 > 
The Nother current and the w;or£; density are thence defined by 

£(L,£,<r) =<F|j fc - 1 i: c a>-z c L 
W(L,£,a) = - < F\ £^a > 

and they satisfy 

d£(L,^a) = W(L,^,a) 

In any natural theory, both the Nother current £(L,£) and the work density 
W(L, £) can be recasted, by suitable covariant integration by parts, as: 

£(L, e, a) = £(L, C, a) + d(W(L, £, <r)) 

(5.19) 

W(L,t,*) = B(L,t,*) + d(£(L,t,*)) 



(5.15) 
(5.16) 

(5.17) 
(5.18) 
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One can easily prove that B(L, £, a) = off-shell (i.e. on any configuration not 
necessarily a solution of field equations), that £(L, £, a) = £(L, £, a) and that 

£, a) vanishes on-shell (i.e. along solutions of field equations). The identity 
B(L,£, a) = is called the generalized Bianchi identity; the current £(L,£,cr) 
is called thereduced current, while U(L,£,o~) is called the superpotential. We 
stress once again that all the objects introduced are algorithmically constructed 
out of the Lagrangian. 

We remark (see (S.18) ) that on shell the Nother current £(L, £, a) obeys 
a continuity equation d£(L,£,cr) = 0, i.e. it is covariantly conserved. The 
conservation laws of this kind are called weak conservation laws since they hold 
on solutions of field equations only, in opposition to the strong conservation 
laws which hold also off-shell. An example of strongly conserved quantity is 
£(L, £, a) — £(L, £, a) which is a closed form along any configuration a , even 
if it is not a solution of field equations (see (5.19) ). 

An example of second order natural theory is General Relativity. Some 
details about this case are found in Section 5 . Further details can be found in 
[2], [28], [27]. The Nother currents, as long as the superpotentials, are used 
in Section 5 to define covariantly conserved quantities in the particular case of 
General Relativity. 

11. Appendix C: Examples 

We shall here discuss some simple examples to illustrate and clarify some 
of the topics introduced above in the paper. 

Let us first consider the conserved quantities of the Schwarzschild solution 
given in its standard form: 

2M\ , o 1 2 . 2 I j/i2 , -„2/)ji2> 



g = -[l ) dt z + -dp 2 + p A (d6H + sin^ 6 d(j) z ) (C.l) 

' ' 2 2M 



P 



For computational convenience we rewrite it in isotropic coordinates (t, r, 9, 4>) : 

s = - jfrw dt2 + 1 1 + Tr ) 4 ( dr2 + r2 (ds2 + sin2 e d * 2) ) (a2) 



where, see [6]: 



2r = -M + p+y/p(p-2M) ^ p= {M + J r) 



2 



(C.3) 
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Let us choose as a background the Minkowski spacetime in isotropic coordinates 



9 = ~ 



(2R-M) 2 
(2R + M) 2 




The background g is in fact a flat metric as one can easily check by direct 
computation of the Riemann tensor. Furthermore, the two metrics g and g 
are matched on the hypersurface B defined by the equation r = R . 

Let us then consider a family of foliations of the region r < R generated 
by the infinitesimal generator 



where e is a (small) parameter. The flow parameter is denoted by s so that, 
if E is a spacelike hypersurface, it can be dragged along the flow <I> S of £ to 
obtain the foliation: 



Since s denotes the affine parameter along the flow of £ , it can be interpreted 
as the time associated to the ADM foliation. We remark that we are in the 
hypothesis of orthogonal boundaries as required throughout the paper. Notice 
also that for e = we recover the ordinary asymptotic time translation £ = d t 
and the ordinary ADM foliation by the hypersurfaces t = constant (see Fig. 2 ) . 



C= (l + e(sin(0) + l))«9* 



(C.5) 



= { s constant} 



s = t/(l + e (sin(0) + 1)) 




e = 



e ^ 



Fig. 2 A 



Fig. 2B 



Let us then choose the ( 3 -parameter) vector field 




(C.6) 
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(where a , (3 and 7 are three real constants) which is a generator of symmetries 
for the first order Lagrangian (1.1) . We remark that £ is a well defined vector 
field on B ; in particular it extends to 9 = and 9 = w . For a = 1 and 
/? = 7 = , £ reduces to the ordinary time translation £ = dt so that we 
expect the corresponding conserved quantity Q£°*[£] to be interpreted as the 
mass of g relative to g on the leaf in the region r < R . We also stress 
that for 7 = the vector field £ is a Killing vector both for g and g . On the 
contrary, for 7 7^ , £ is not Killing both for g and g . In the case 7 7^ , 
£ has not a direct physical interpretation though it is a symmetry generator 
for the Lagrangian (1.1) and it algorithmically generates the covariant Nother 
conserved quantity Qi^tC] • We shall use it to illustrate how time conservation 
may be related to covariant conservation along different foliations. 

If we calculate the Nother conserved quantity (according to eq. (6.1) with 
2k = 16-7T , i.e. in geometric units G = c = 1 ) we get the following result 



n \t 1 M n nf-t M V 
Qf£ e Ki Q\ = a 7 — a 1 s e 

^ r , , / M ttM \ / M\ 
Q E! [^]=^ T -7 1 6-"j 

nTtm /„, M 2 \ /ttM 2 \ 



Thus, when £ is the Killing vector d t (a = 1 , j3 = 7 = ) the relative 
mass 

Ql?[d t ] = M - ^ (C.8) 

is time-conserved along any foliation of the family generated by £ , since it 
does not depend on the affine parameter s . We also stress that the conserved 
quantity (i.e. letting R tend to infinity) is always Q™[dt] = M , i.e. it reduces 
to the expected value for total mass. 

If 7 7^ the vector £ is not Killing. Despite of the fact that the conserved 
quantity is not time-conserved in general, it is still time-conserved along a 
particular foliation (namely, the one corresponding to e = , i.e. ( = dt ; see 
Fig. 2 A). We remark that in this case the flow of the Nother current £[£] = 
£^ ds^ through <9E S vanishes even if £ M n M | B 7^0. 

On the contrary, when e 7^ (see Fig. 2 B) the Nother charge explic- 
itly depends on the time s , i.e. on the particular leaf on which it is computed. 
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We finally remark that in any case the quantity <5e°*[C] given by equation 
(C.7) is covariantly conserved (as any Nother conserved quantity is) and it 
consequently obeys a continuity equation. 

To end up this first example let us now compute the quasilocal energy 
(6.15) . We consider the case depicted in Fig. 2 A, i.e. the foliation induced by 
( = d t . We choose, as infinitesimal symmetry £ , the unit timelike normal u : 



2r + M 
U= ^~M dt 



(C.9) 



If we calculate the Nother conserved quantity (6.1) we now obtain the following 
result: 

1 M(2R + M) 

Qv t [u,g] = j 2R _ M 



<5s t [u, g] 

qi°;[u] = m + 



1 (2R 2 - MR - M 2 )M 



2 R(2R - M) 
1 M{4R 2 - M 2 ) 
4 (2R-M) 2 
M 2 



(CIO) 



2R 



In spherical coordinates (t,p,0,4>) the fourth expression of (CIO) may be 
rewritten as: 



QTM = po 



l 



l- 



2M 
Po 



(Cll) 



where 2R = —M + po + \J po(po — 2M) ; see (C.3) . As expected, expression 
(C.ll) perfectly agrees with the value of the energy computed in [3], formula 
(6.14) . 

Let us now consider another example of conserved quantity in a finite region 
D . It is a completely different example since it does not require the match of the 
solutions on the boundary B of the finite region under consideration. Neither 
the condition of orthogonal boundaries is here required. 

Let us consider the Kerr spacetime in ingoing Kerr-Schild coordinates 
(*, r, 6, 4>) , given by: 



g = g + 2Mrp~ 



dt + dr — a sin 6>d 4> 



(C.12) 
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where p 2 = r 2 + a 2 cos 2 9 , M 2 > a 2 . Let us choose the flat background g as: 



g = — dt 2 + dr — a sin 



2 9d(P +p 2 d9 2 + sin 2 9 d 2 



(C.13) 



The metrics <7 and g are matched at infinity. Let us consider the regions D 
inside the hypersurface B defined by r = R and the ADM foliation E t = 
{t constant} generated by the vector field d t . We stress that g and # do not 
match on B unless R is let to tend to infinity. 

Let us finally choose as symmetry generator the ( 2 -parameter) vector 



which is a Killing vector for g and g ( a and (3 are two real constants) . 
The Nother conserved quantities one obtains are 



which reproduce the expected values of the relative mass and of the angular 
momentum in the region t = constant and r < R . Notice that the result 
is independent on R meaning that all the energy and angular momentun is 
"buried in the singularity". We remark that setting anywhere a = the 
Schwarzschild solution is recovered. The relative mass M we obtain in this 
case (see equation (C.15) with a = 1 , (3 = ) does not agree with the value 
found above (see equation (C.8) ) because of the two different matches selected. 

We remark that the quantities (C.15) are also time-conserved (even if the 
metrics are not matched at B). In fact the flow integral j B S vanishes since 
£ is tangent to B and it is a Killing vector of both g and g (see Section 7). 
Notice also that £ a n a = , i.e. the Nother current has no flow through any 
part of B . Consequently, the associated conserved quantity is time conserved 
along any ADM foliation of the region D . 

[The calculations in this Appendix have been carried out by using tensor 
package of MapleV, see [43]. They are the very direct application of formula 
(6.1) , just computed on the configuration (g,g) .] 



£ = a d t + (3 



(C.14) 



Ql ot [£] =aM -/3 Ma 



(C.15) 
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